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Abstract. We have used the Maranol field observations 
with ISOPHOT at 170 /xm to search for the Cosmic Far- 
InfraRed Background fluctuations. This field is part of the 
FIRBACK project (Puget et al. 1999). For the first time, 
fluctuations due to unresolved extra-Galactic sources are 
isolated. The emission of such sources clearly dominates, 
at arcminute scales, the background fluctuations in the 
lowest Galactic emission regions. 

The study presented here is based on a power spectrum 
analysis which allow to statistically separate each of the 
background components. With this analysis, we clearly 
show that wc detect the CFIRB fluctuations. 

Key words: Cosmology: diffuse radiation - in- 
frared:ISM:contimium - infrared: general 



1. Introduction 

The major component of the extra-Galactic background 
is formed by the integrated light of all distant galaxies 
that are not resolved in the beam of the instrument. The 
detection of the far-infrared part of this background (the 
CFIRB, Cosmic Far-InfraRed Background) with COBE 
(Puget et al. 1996; Fixsen et al. 1998; Hauser et al. 
1998; Lagache et al. 1999) has opened new perspectives 
in our understanding of galaxy evolution. The energy 
contained in the CFIRB, compared to the one in the 
UV/optical/near-IR domain (Lagache 1998; Gispert et al. 
in preparation) shows that a new population of galaxies 
has to be considered in the general framework of galaxy 
evolution: the far-infrared galaxies. The closest and/or 
most luminous ones have been discovered by IRAS (see 
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for a review Sanders & Mirabel 1996). Nowadays, the IR 
galaxies are hunted down by several deep cosmological sur- 
veys such as the ones with SCUBA at 850 /im (Hughes et 
al. 1998; Small et al. 1997, 1998; Lilly et al. 1999; Barger 
et al. 1999), ISOPHOT at 170 /xm (Kawara et al. 1998; 
Puget et al. 1999) and ISOCAM at 15 ^^m (Oliver et al. 
1997; Aussel et al. 1999; Desert et al. 1999; Elbaz et al. in 
preparation) . 

FIRBACK is a deep cosmological survey dedicated to 
the investigation of the nature of the galaxies that con- 
tribute to the CFIRB at 170 ^m (Dole et al. 1999) with 
the ISOPHOT instrument (Lemke et al. 1996). It covers 
about 4 square degrees. The first FIRBACK field to be 
observed and analysed is the Maranol region centered on 
RA(J2000)=3hl3m9.6s and DEC(J2000)=-55d03m43.9s. 
This field covers 0.25 square degree and contains 24 extra- 
Galactic sources with fluxes ranging from 100 to 850 mJy 
(Puget et al. 1999). It lies in a very low cirrus emis- 
sion region with HI column density equal on average 
to 10^° H cm^^. This low cirrus contamination allows 
to search for brightness fluctuations due to unresolved 
sources. 

In this paper, we use the FIRBACK Maranol field to 
isolate the background fluctuations with a high signal 
to noise ratio. We then discuss the separation of the 
background fluctuations between the Galactic and extra- 
Galactic component. 

Throughout this paper, "background" refers to the extra- 
Galactic and Galactic backgrounds, and "CFIRB" to the 
extra-Galactic background at far-IR wavelengths. 



2. Sources of the CFIRB 

The CFIRB is made of sources with number counts as a 
function of flux which can be represented, for the present 
discussion, by a simple power law: 



iV(> S) = No 



(1) 
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For an Euclidian Universe, a=1.5. Obviously, these num- 
ber counts need to flatten at low fluxes to insure a finite 
value of the background. Thus, we assume in the present 
discussion that a—0 for S < S* . 

The intensity of the CFIRB, induced by all sources 
with flux up to Smax, IS givcn by: 



iCFIRB 



dN 
'dS 



SdS 



(2) 



For the simple Euclidian case {a— 1.5), the CFIRB inte- 
gral is dominated by sources near S*. 

Fluctuations from sources below the detection limit 
So (which corresponds to the flux of sources either at the 
confusion or at the sensitivity limit) are given by: 

^2^ ['" S'^dS JyVsr (3) 

Jo "'-' 

Computing ^ using Eq. (1), this gives: 



-A^o^o' 



So 



2-a 



JyVsr 



(4) 



For the simple Euclidian case , the CFIRB fluctuations are 
dominated by sources which are just below the detection 
limit So- 

Nevertheless, it is well known that strong cosmological 
evolution, associated with a strong negative K-correction, 
could lead to a very steep number count distribution (see 
for example Guiderdoni et al. 1998; Franceschini et al. 
1998). For a >2, the CFIRB integral is still dominated 
by sources near S* but its fluctuations are now also dom- 
inated by sources close to S*. In such a case, observations 
of the fluctuations will constrain sources that dominate 
the CFIRB. Futhermore, investigations of problems like 
the fraction of AGN contributing to the CFIRB, which 
has been identified as a major question of the coming 
years, cannot be conducted only with individually de- 
tected sources. These observed sources correspond to the 
bright part of the luminosity function and only represent a 
few percent of the CFIRB. Moreover, the fraction of AGNs 
could very well be increasing with luminosity (Sanders & 
Mirabel 1996). Thus the only way to tackle this prob- 
lem will be to study the correlation between the extra- 
Galactic background fluctuations in the X-rays and Far- 
IR domains. 




Fig. 1. Continuous line; ISOPHOT Saturn footprint power 
spectrum (Ws^). Dotted line: ISOPHOT model footprint 
power spectrum (Wnifc). The difference in solid angle between 
the two footprints is equal to 25%. 



3. Summary on the data reduction and calibration 

The data reduction and calibration for this field have been 
detailed in Lagache (1998), Puget et al. (1999) and Puget 
& Lagache (1999). We are just going to summarise here 
the different steps. 

We use PIA, the ISOPHOT interactive analysis soft- 
ware version 6.4 (Gabriel et al. 1997), to correct for instru- 
mental effects, glitches induced by cosmic particles and 
to provide an initial calibration. First we apply the non- 
linearity correction. Deglitching is performed for each in- 
dividual ramp and then the mean signal per position is de- 
rived by averaging the linear fit on each ramp (see Gabriel 
et al. 1997 for details). The signal shows a long term drift 
which is corrected by using the two internal calibration 
(FCS) measurements bracketing each individual observa- 
tion. Inside each raster, the long term drift represents less 
than 10% of the signal. The calibration is performed by 
deriving the mean value of the two FCS measurements. 
The contribution of the long term drift is thus lower than 
5%. Flat fielding is achieved using the high redundancy in 
the measurements. Glitches inducing long term drifts are 
also corrected. The fiux is finally projected on a 10"xlO" 
coordinate grid using our own projection procedure. 



In the Far-IR, present observations show a very steep 
slope a=2.2 (Dole et al. 1999). Sources detected above 
So contribute only to less than 10% of the CFIRB (with 
So=120 mjy). Thus, it is essential to study the extra- 
Galactic background fluctuations which are likely to be 
dominated by sources with a flux comparable to those 
dominating the CFIRB intensity. 



The final calibration is achieved using the footprint 
measured on Saturn. It is presented in Fig. 1 of Puget et 
al. (1999). Brightnesses derived using this footprint are 
in very good agreement with DIRBE. The footprint mea- 
sured on Saturn has been compared to the one derived 
from the diffraction model, including the tripode (Klaas 
et al. in preparation). The model takes into account the 



G. Lagache et al.: Detection of the extra-Galactic background fluctuations at 170 /im 



MJy/s 



MJy/sr 




13 12 

a2nnn (h^) 



Fig. 2. a. The Maranol region observed with ISOPHOT at 170 fim. b. same as (a) with detected sources removed. 



bandpass filter and is computed for a i^Iiy=constant spec- 
trum. Thus, comparison between the model and Saturn 
footprints allows to test the wavelength dependence of 
the footprint. The solid angle supported by the two foot- 
print differs only by 25%. Comparison of the two foot- 
prints power spectra, Wsj^ (Saturn) and Wmj^ (model), 
is shown on Fig. 1. The Saturn FWHM footprint is quite 
larger than the model one. The difference between these 
power spectra will be considered in Sect. 5.2. 

The Maranol observations consist of four 19x19 
rasters. Each raster was performed in the spacecraft (Y,Z) 
coordinate system which is parallel to the edges of the de- 
tector array. The field area covered by each raster is about 
30'x30'. The exposure time is 16s per pixel. Rasters were 
performed with one pixel overlap in each Y and Z direc- 
tion. Therefore, for the maximum redundancy region, the 
effective exposure time is 256 sec per sky position (16 sec 
per pixel, 4 rasters and redundancy of 4). Displacements 
between the four raster centers correspond to 2 pixels. 
This mode of observation, which does not provide proper 
sampling of the point-spread-function, was chosen delib- 
erately because it allows a very clean determination of the 
instrumental noise. 

The final map is made of 16 coadded independent 
maps. The extraction of extra-Galactic sources and de- 
termination of their fluxes are made on this final coadded 
map. To compute the instrumental noise, we have mea- 
sured the flux of the detected sources on the 16 indepen- 
dent maps. We then compute for each source the stan- 
dard deviation with respect to the flux measured on the 



coadded map. The average of all the standard deviations 
gives the instrumental noise. We obtain 1.74 mJy rms in 
a 89.4" X 89. 4" pixel. The signal in the map is around 
3 MJy/sr. This leads to a very high signal to noise ra- 
tio in the map of about 300. 

Using the four maps of the Maranol field presented 
in Fig. 2 of Puget et al. (1999), we can make an esti- 
mate of the instrumental noise power spectrum, making 
the difference between the independent maps. We know 
that this power spectrum contains low spatial frequency 
noise, due to small distorsions in the flat field and slow 
response changes during the observations. With such map 
differences, we obtain instrumental noise power spectra 
for pairs of maps which are all more than 10 times smaller 
than the measured total flux power spectum (which is pre- 
sented in Fig. 3). The instrumental noise power spectrum 
is rather flat in the range k=:[0.1, 0.5] arcmin"^ (about 
30 Jy2/sr at k--0.2 arcmin^^ and 10 Jy^/sr at k-0.4 
arcmin^^) and is in very good agreement with the one de- 
duced from the noise (assumed to be white) computed on 
the extra-Galactic sources (1.74 mJy rms in a 89.4" x89.4" 
pixel). The measured instrumental noise power spectrum 
is subtracted to the measured power spectrum before the 
analysis. 

4. The extra-Galactic sources 

In the Maranol field, 24 extra-Galactic sources are de- 
tected with fluxes ranging from 100 to 850 mJy (Puget 
et al. 1999). Sources have been extracted up to the con- 
fusion limit (3(7=67.2 mJy) and source fluxes been com- 



G. Lagache et al.: Detection of the extra-Galactic background fluctuations at 170 /im 



puted using the Saturn footprint measurement for which 
a significant fraction of the integral is in the wing (and 
which is in good agreement with the footprint model). 
The source subtraction from the Maranol map is very 
straightforward. We subtract at each source position the 
footprint, sampled at 10 arcsec and normalised to unit in- 
tegral, multiplied by the source flux. The map obtained 
after the source subtraction is given Fig. 2 together with 
the original map. The standard deviation in the source 
subtracted map is around 0.14 MJy/sr, which is very close 
to the one of the original map which is about 0.15 MJy/sr. 
For the mean brightness, the difference between the two 
maps is less than 1%. These numbers clearly show that 
the background fluctuations are not dominated by strong 
sources detected above the confusion limit. As for the orig- 
inal map, the source subtracted map has a very high signal 
to noise ratio. 

Background fluctuations (Fig. 2b) are made of two 
components: the Galactic cirrus one and the extra- 
Galactic one (the instrumental noise is negligible and is 
not correlated). We are going now to discuss the relative 
contribution of these two components. 



5. Extra-Galactic and Galactic fluctuation separa- 
tion 

Our separation of the extra-Galactic and Galactic fluctu- 
ations is based on a power spectrum decomposition. This 
method allows us to discriminate the two components us- 
ing the statistical properties of their spatial behaviour. 

We know from previous works that the cirrus far- 
infrared emission power spectrum has a steep slope in k^^ 
(Gautier et al. 1992; Kogut et al. 1996; Herbstmeier et al. 
1998; Wright 1998). These observations cover the relevant 
spatial frequency range. Futhermore, Abergel et al. (1999) 
have shown that this spatial structure, measured from ex- 
tinction, extend to arcsecond scales. Despite a large obser- 
vationnal effort, no characteristic scale has been identified 
in the diffuse interstellar medium (Falgarone 1998). 
For the extra-Galactic fluctuations, little is known. The 
power spectrum of galaxy clustering has been mainly stud- 
ied from optical surveys. The main result is that the power 
spectrum can be represented by a simple power law with 
an index equal to -1.3 at spatial frequency larger than ~0.1 
degree" 1 (Groth & Peebles 1977; Maddox et al. 1990; Pea- 
cock 1991; Kashlinsky 1992). For small spatial frequency, 
if we assume the standard inflationary model (f2=l and 
the Harrison- Zeldovich power spectrum on large scales), 
the galaxy power spectrum is proportional to k. At Far-IR 
wavelengths, nothing is known about the galaxy cluster- 
ing. A spatial and redshift distribution similar to the opti- 
cal surveys (as the "Automatic Plate Measuring" survey) 
will give a negative slope around -1.3 whereas a Poisso- 



nian distribution will give a flat power spectrum. 

In summary, we know that the cirrus emission power 
spectrum is steep and does not present any characteris- 
tic scale. The extra-Galactic component is unknown but 
certainly much flatter. We thus conclude that the steep 
spectrum observed in our data at k<0.15 arcmin"^ (Fig- 
3) can only be due to the cirrus emission. The break in 
the power spectrum at k~ 0.2 arcmin~^ is very unlikely 
to be due to the cirrus emission itself which is known not 
to exhibit any prefered scale. Thus, the power spectrum 
is a powerful tool for a flrst order separation between the 
Galactic and extra-Galactic background component. We 
now give and discuss the quantative separation of the two 
components. 

5.1. Normalisation of the cirrus emission fluctuations 

The cirrus emission has a filamentary structure which, as 
discussed above, exhibits a spatial power spectrum which 
has been first determined at 100 fim from IRAS data by 
Gautier et al. (1992): 



Pc. 



m 



1.4 X IO-i^BqS I _ J jy2/g^ 



-3 



(5) 



with kg = 10"^ arcmin"^, Bq the cirrus brightness at 
100 /J.m in Jy/sr and k the spatial frequency in arcmin"^. 
The main problem in this determination comes from the 
normalisation which has been established in high cirrus 
emission regions. This normalisation does not necessarily 
apply to our very low cirrus emission field. For this reason, 
the normalisation is directly determined using the mea- 
sured power spectrum (Fig. 3). Using the low frequency 
data points and assuming a k^^ dependance, we find a 
power spectum of: 



^cc,i7o(fc) = 7.9+^:^ k-' JyVsr 



(6) 



with k in arcmin"^. The main uncertainty comes from the 
footprint used for deriving -Pcc,i7o(fc)- The normalisation 
of the power spectrum derived from the low frequency data 
points agrees reasonably well with Wright (1998) DIRBE 
results at high Galactic latitude. It is also in good agree- 
ment with the Gautier et al. (1992) normalisation scaled 
to 170 /im (Appendix A). 

5.2. Measured power spectrum 

We can now compare the power spectrum measured on 
the source subtracted map (Fig. lb) with our estimated 
one for the non cosmological components defined as: 



Poc 



est 



Pccayo X Wy. 



(7) 



where Wy^ is the ISOPHOT 170 /im beam power spec- 
trum, and Pec, 170 the cirrus confusion noise. Fig. 3 shows 



U) 



>^ 



CL 



10 



5 







4 







3 







2 







G. Lagache et al.: Detection of the extra-Galactic background fluctuations at 170 /im 



r '^-^ 



^ 



'■. ^ 



\ 
\ 



>^ )^ 



)K 



*; 



,^ 



>K^^)K 



>K 



)l^ 



0.1 



.0 



k (arcmin ) 



Fig. 3. Power spectrum of the source subtracted map (*). The dashed and dotted lines correspond to the cirrus confusion noise 
multiplied by the footprint power spectrum (Wsu and Wmi^ respectively) 



the measured power spectrum together with Pggi; com- 
puted for W]j=Wsjj. and W]j=Wm]j (as defined in Fig. 1). 
We clearly see an excess between k=0.25 and 0.6 arcmin"^ 
which is more than a factor 10 at k=0.4 arcmin"'^. Any 
reasonable power law spectrum for the cirrus component 
multiplied by the footprint leads, as can be easily seen 
in Fig. 3, to a very steep spectrum at spatial frequency 
k>0.2 arcmin^ ^. This is very different from the observed 
spectrum. The excess is quite independent of the footprint 
used. Moreover, it is more than 10 times larger than the 
measured instrumental noise power spectrum. Therefore, 
as no other major source of fluctuations is expected at this 
wavelength, the large excess observed between k=0.25 and 
0.6 arcmin^^ is interpreted as due to unresolved extra- 
Galactic sources. The anisotropics of the background are 
thus dominated, at angular scales close to the angular res- 
olution, by the structure of the CFIRB. This is the first 
detection of the fluctuations of the CFIRB. 



dynamic range to constrain the slope of the CFIRB fluc- 
tuation power spectrum: our present study docs not allow 
to clearly constrain the clustering of galaxies. However, 
the power spectrum mean level can be determined. 
Power spectra with simple power law with index smaller 
than -1 arc not compatible with our measured power spec- 
trum. For indexes in the range [-1,0], the different power 
laws give equivalent results. For example, on Fig. 4 is 
shown the total power spectrum obtained by summing 
up the cirrus confusion noise and a constant CFIRB fluc- 
tuation power spectrum. Pes ~ 7400 Jy^/sr. This white 
noise power spectrum level is in very good agreement with 
the one predicted by Guiderdoni et al. (1997). This gives 
CFIRB rms fluctuations around 0.07 MJy/sr (for a range 
of spatial frequency up to 5 arcmin^^). These fluctuations 
are at the ^9 percent level, which is very close to the pre- 
dictions of Haiman & Knox (1999). 



CFIRB fluctuations are clearly detected between 
k=0.25 and 0.6 arcmin^^. This interval has not enough 
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Fig. 4. *: power spectrum measured on the map. Dotted line: 
total power spectra obtained by summing up the cirrus confu- 
sion noise (adjusted to our large scale points) and a flat CFIRB 
fluctuation power spectrum Pes ~ 7400 Jy^/sr (we have used 
here the Saturn footprint power spectrum). 



6. Conclusion 

We have shown in the Maranol region (30'x30') that the 
extra-Galactic background fluctuations are well above the 
instrumental noise and the cirrus confusion noise. The ob- 
served power spectrum shows a flattening at high spa- 
tial frequencies which is due to extra-Galactic unresolved 
sources. A preliminary study on the FIRBACK Nl field 
(2 square degrees) show exactly the same behaviour. The 
statistical analysis of the background fluctuations using 
the full FIRBACK fields, covering different high Galactic 
latitude regions, will allow to constrain models of galaxy 
evolution and will be a very useful tool for cosmological 
studies until large passively cooled telescopes in space be- 
come available which will allow number counts down to a 
few mjy. The next step consists of removing the cirrus con- 
tribution using gas tracers (the Hq and the 21cm emission 
lines) to isolate the extra-Galactic fluctuation brightness 
and thus constrain the clustering of galaxies. 
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Appendix A: Comparaison of the measured cir- 
rus power spectrum with the Gautier et al. (1992) 
normalisation: 



The power spectrum of the cirrus emission derived by 
Gautier et al. (1992) from IRAS data at 100 fim follows: 



Pcc,ioo(fc) = 1.4 X 10-^^Bo' ( f- ) Jy'/sr 



(8) 



with kg = 10 ^ arcmin ^, Bq the cirrus brightness at 
100 /ini in Jy/sr and k the spatial frequency in arcmin^ ^. 

The far-infrared Galactic sky at high latitudes has es- 
sentially three components: the zodiacal background, the 
emission from Galactic cirrus and the emission from in- 
frared galaxies at all redshifts. The average total flux 
in our map is 3.16 MJy/sr. The zodiacal emission has 
been derived at the time of our observation using the 
model of Reach et al. 1995 (calibrated on DIRBE data) 
and is equal to 0.67±0.07 MJy/sr. Following Lagache et 
al. (1999), the CFIRB at 170 fim is equal to 0.77±0.25 
MJy/sr. Thus, we obtain for the cirrus emission at 
170 ^m, I(170) = 1.72±0.26 MJy/sr. Using the color ratio 
for the HI diffuse emission I(170)/I(100)=1.96 (Lagache 
et al. 1999), we derive Bo= 0.88±0.13 MJy/sr= 8.8±1.3 
10'"' Jy/sr and the normalisation of the power spectrum at 
170 fim becomes: 



1.4 X 10" 



/(170) 



1 2 



/(lOO) 



B^ ^ 3.63+^^ 



Using Eq. (8), this gives: 

PcCA70{k) = Smtli X fc-3 Jy2/g^ 

with k in arcmin^ ^. 



X 10^ JyVsr (9) 



(10) 



In fact Gautier et al. results flt to within a factor of 
2 the Eq. (8). Taking into account this uncertainty, the 
cirrus power spectrum becomes: 



Pcc,i7o(fc) = 3.63+^:t? X k-" JyVsr 



(11) 



We do not expect in such a low column density field 
(Nhj ^ 10^" H cm^^) any contribution of dust associ- 
ated with molecular gas (Lagache et al. 1998). Therefore, 
we can argue that the normalisation of the cirrus confu- 
sion noise at 170 /ini, for our field, is between 1.1 and 11.1 
Jy^/sr. 

Our normalisation computed from the low frequency 
data points is around 7.9 Jy^/sr which is in the range 
of Gautier et al. results. This value of 7.9 Jy^/sr agree 
also reasonably well with Wright (1998) measures using 
DIRBE data. 
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